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Abstract 

We developed the generalized tight-binding model to study the magneto-electronic 
properties of AAB-stacked trilayer graphene. Three groups of Landau levels (LLs) are 
characterized by the dominating subenvelope function on distinct sublattices. Each 
LL group could be further divided into two sub-groups in which the wavefunctions are, 
respectively, localized at 2/6 (5/6) and 4/6 (1/6) of the total length of the enlarged 
unit cell. The unoccupied conduction and the occupied valence LLs in each sub¬ 
group behave similarly. For the first group, there exist certain important differences 
between the two sub-groups, including the LL energy spacings, quantum numbers, 
spatial distributions of the LL wavefunctions, and the field-dependent energy spectra. 
The LL crossings and anticrossings occur frequently in each sub-group during the 
variation of field strengths, which thus leads to the very complex energy spectra and 
the seriously distorted wavefunctions. Also, the density of states (DOS) exhibits 
rich symmetric peak structures. The predicted results could be directly examined by 
experimental measurements. The magnetic quantization is quite different among the 
AAB-, AAA-, ABA-, and ABC-stacked configurations. 
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1. Introduction 


Graphene, being a one-atom-thick layer of carbon atoms densely packed in a two- 
dimensional honeycomb lattice, has attracted a lot of theoretical and experimental re¬ 
search [MI]. It exhibits many unusual physical properties, e.g., a rich magnetic quan¬ 
tization M, half-integer Hall effect [5HT5], high Young’s modulus HMD, high Fermi 
velocity (10 6 m/s), and others. Graphene could play an important role in technological 
applications such as electric circuits [221123] . field-effect transistors [24j,[25|, light-emitting 
diodes [251 [281 . solar cells [25H32], and durable touch screens [551154] . Due to the hexagonal 
symmetry with a rotational angle of 60°, a non-doped graphene is a zero-gap semiconduc¬ 
tor with a vanishing density of state at the Fermi level. The essential electronic properties 
can be drastically changed by the layer number [3-51 137j . stacking configuration [5TH42] . 
magnetic held |431i45], electric held [45114?] . dopping j48[[49j, mechanical strain [5Q1T52] . 
and temperature variation [531'-TT . Few- and multi-layer graphenes have been success¬ 
fully produced by experimental methods such as exfoliation of highly orientated pyrolytic 
graphite [55H58] . metalorganic chemical vapour deposition (MOCVD) [5TH55] . chemical 
and electrochemical reduction of graphene oxide [5TH55] . and arc discharge mm There 
exist important stacking configurations, including AAB [57] [58], [69], ABC [ 59H60II66H69] . 
AAA [521153] , ABA Qp.5j -JiDlTin. and twisted [52)[5Sj and turbostratic ones [54] . The in¬ 
terlayer atomic interactions and stacking configurations induce the rich electronic properties 
of graphene. In this work, we investigate the complex relationship between the magnetic 
quantization and the interlayer atomic interactions in AAB-stacked trilayer graphene by 
using the generalized tight-binding model. 

The low-lying energy dispersions of monolayer graphene possess a pair of linear bands 
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intersecting at E F . These energy bands are further quantized by the application a uni¬ 
form perpendicular magnetic field B = Bqz [72]. The magneto-electronic spectrum of the 
isotropic Dirac cones satisfies a simple E c,v oc Ey/nB^ relationship, where n is the quantum 
number, and c and v denote the conduction and valence Landau levels (LLs), respectively. 
The band structure and the LLs of monolayer graphene have been verified by a number of 
experimental methods [531I751I7I] . Bilayer graphene, being held together by Van der Waals 
interactions, can exhibit the highly symmetric AA and AB configurations. The former and 
the latter have, respectively, two pairs of linear and parabolic bands. The two isotropic 
Dirac-cones in the AA-stacked system are magnetically quantized into two groups of LLs 
with monolayer-like wavefunctions and energy spectrum m Also, the AB-stacked sys¬ 
tem presents two groups of LLs, with each LL having a single-mode wavefunction in the 
absence of LL anticrossing [77]. The field-dependent LL energy spectrum of AB-stacked 
bilayer graphene is different from that of monolayer graphene. It includes a few intergroup 
LL anticrossings at large field strengths (B 0 > 100T). The band structures [7511751175] and 
the first LL groups have been experimentally verified for bilayer graphene. 

The magneto-electronic properties of AAA- and ABA- stacked trilayer graphenes can be 
regarded as combinations of those of monolayer graphene [35], and AB-stacked bilayer and 
monolayer graphenes [383, respectively. Distinctly, the zero-field band structure of ABC- 
stacked trilayer graphene exhibits three pairs of energy dispersions: linear cones, sombrero- 
shaped bands, and parabolic bands 0138]. Accordingly, the I?o-dependent energy spectrum 
has both intergroup and intragroup LL anticrossings, in which the latter are arised from 
the sombrero-shaped energy bands. The band structures and the first LL groups of these 
stacking systems have been examined by numerous experimental measurements [681179] 18T] . 
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As for the stacking symmetry, the AAB-stacked system is in sharp contrast with the AAA- 
and ABA-stacked ones, but is also different from the BAC-stacked one. Clearly, the AAB- 
and ABC-stacked trilayer graphenes can not be regarded as a superposition of two or three 
sub-systems. 

We develop the generalized tight-binding model, based on the subenvelope functions 
on the distinct sublattices, to study the rich electronic properties of AAB-stacked trilayer 
graphene. This work shows that the AAB-stacked configuration presents an abnormal 
band structure, in which there are pairs of oscillatory, sombrero-shaped and parabolic 
bands, which are different from those of the other stacking systems mentioned above. 
The special band structure could be verified by angle-resolved photoemission spectroscopy 
(ARPES) [82]. The LL evolution under a magnetic field reveals a complex pattern of 
LL anticrossings and splittings, as a result of the specific interlayer atomic interactions 
derived from the full tight-binding model. The magnetoelectronic spectra are directly 
reflected in DOS which can be examined by experimental measurements using scanning 
tunneling spectroscopy (STS). The R 0 -dependent energy spectra exhibit intragroup and 
intergroup LL anticrossings. Especially, the low-lying LLs anticross frequently, leading to 
the existences of non-well-behaved and perturbed LLs. The spatially dramatic changes of 
the LL wavefunctions in the anticrossings can be verified by scanning tunneling microscope 
(STM) measurements [571158], as done for 2D electron gas [83] and topological insulators 
[MJ. In this paper, the important differences among trilayer graphenes such as their energy 
band structures, LL splittings, LL ordering, and LL anticrossings, are investigated. 
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2. The generalized tight-binding model 


The low-energy ^-electronic structure of AAB-stacked trilayer graphene, mainly coming 
from the 2 p z orbitals, is calculated with the generalized tight-binding model. The two 
sublattices in the Zth (l = 1,2,3) layer are denoted as A 1 and B l . The first two layers, 
shown in Fig. 1, are arranged in the AA-stacking configuration; that is, all carbon atoms 
have the same (x,y) projections. The third layer can be obtained by shifting the first (or the 
second) layer by a distance of b along the armchair direction. In this system, the A atoms 
(black) have the same (x,y) coordinates, while the B atoms (red) on the third layer are 
projected at the hexagonal centers of the other two layers. The interlayer distance and the 
C-C bond length are, respectively, d = 3.37/1 and b = 1.42/1. There are six carbon atoms in 
a primitive unit cell. The low-energy electronic properties are characterized by the carbon 
2 p z orbitals. The zero-field Hamiltonian, which is built from the six tight-binding functions 
of the 2 p z orbitals, is dominated by the intralayer and the interlayer atomic interactions 
There exist 10 kinds of atom-atom interactions corresponding to the 10 atomic hopping 
integrals which appear in the Hamiltonian matrix. y 0 = —2.569 eV represents the nearest- 
neighbor intralayer atomic interaction; 74 = —0.263 eV, 72 = 0.32 eV present the interlayer 
atomic interactions between the first and second layer; 73 = —0.413 eV, 74 = —0.177 eV, 
75 = —0.319 eV are associated with the interlayer atomic interactions between the second 
and third layer; y 6 = —0.013 eV, 77 = —0.0177 eV, and 73 = —0.0319 eV relate to 
the interlayer atomic interactions between the first and third layer; and 79 = —0.012 eV 
accounts for the difference in the chemical environment of A and B atoms. The hopping 
integrals 71 , 73, and 75 belong to the vertical interlayer atomic interactions, while the others 
are non-vertical ones. 
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When applying a uniform perpendicular magnetic field, the unit cell becomes enlarged 
as indicated in Fig. 1. There appears an extra Peierls phase Gr in the tight-binding func¬ 
tions. Gr has the form of f^A.dl, where A is the vector potential, and </> 0 = hc/e is 
the flux quantum. The vector potential in the Landau gauge is chosen as A = (0,Bx,0). 
The Peierls phase has a period of 20o/0 = 2 Rr along the x-axis. Under the effect of 
the magnetic held, the unit cell becomes enlarged as a rectangle with 12 R r atoms in¬ 
cluded. There are 12 Rr tight-binding functions, which are arranged in the sequence 
{l^i l^i); \Bf), |-4?), | Bf ),..., \A\ Rb \B\ Rb ), 1 ^ 2 ^), |-Bf Rs ), \A\ Rb ), |-B? Rs )}. The su¬ 
perscript Hamiltonian is a 12 R r x 12 R r matrix, in which the non-zero matrix elements 
can be presented by the equations below. 


(B)\H\A\) = 7o ^2^exp[ik.(R^ 



2iir 

4>o 




— 7o + yoOWy+i 


= 7o ti,i<7ij + loqcri, j+ i 

(Bj\H\Af) = 7 0 t 3 ) i^ij-i + 7o q<Ti,j 

(A)\H\A]) = (B]\H\B]) = 


( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 


(B)\H\A$) = (A)\H\B*) = 


( 6 ) 
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{ A ]\ H \ A V ) = 73^ ( 7 ) 

(Sj|fT|S?)=7fet 2ii (7y (8) 

(5JI5IA) = 'yrfijCTij + 7 4 g<7ij+i (9) 

+ 74?^ij (10) 

(^1^14) = 76<Ty (11) 

(S?|i/'|S^> = 78^2, i<Tij + 78^7+1 (12) 

(5)1/71 A) = 77^i,i^j + 7r9^*j+i (13) 

( A )\ H \ B i) = + 77 q(Tij (14) 


The four independent phase terms are: 
tij = exp{i[-(k x b/2) - (y/3k y b/2) + 7r</>(j - 1 + 1/6)]} 
+ exp{i[-(k x b/2 ) + (• \/2>k y b/2 ) - %(f)(j - 1 + 1/6)]} 
t 2 ,j = exp{i[-(k x b/ 2) - (y ^3k y b/2) + 7r0(j - 1 + 3/6)]} 
+ exp{i[— (k x b/2) + (\/3 k y b/2) — ir(j)(j — 1 + 3/6)]} 
fij = exp{f[-(/c a; fe/2) - (V3kyb/2 ) + 7T0(j - 1 + 5/6)]} 
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+ exp{i[-(k x b/ 2) + (\/3k y b/2) - n(j)(j - 1 + 5/6)]} 
q = exp{ik x b}. 

In order to deal with the small values of magnetic field strength and a huge Rb, we arrange 
the Hamiltonian in a band-like symmetric matrix. The Landau wave functions, which are 
investigated to identify the spatial distributions of the Landau levels, can be expressed as 

3 2 Rg 

W = \A‘J + B L\B' m )), (15) 

1=1 m— 1 

where A l m (B l m ) are the subenvelope functions presenting the amplitude of the tight-binding 
functions based on the mtli A (B) atom at the Zth layer in the unit cell. 

3. The zero-field band structure and quantized Landau levels 

The zero-field band structure of AAB-stacked trilayer graphene consists of three pairs of 
conduction and valence subbands labeled S[’ v } S^, and S^’ v , as shown by the solid curves 
in Fig. 2(a). Near the Fermi energy, the two subbands which belong to the first pair, Sl’ v , 
have strong oscillatory energy dispersions. The conduction subband starts to increase from 
a local minimum value of about 4 meV at the K point (the corner of the first Brillouin 
zone), along the KM and Kr directions. After reaching a local maximum value of about 58 
meV, it decreases until reaching a local minimum energy again (4 meV), and then grows 
steadily. The curvature of the valence subband is in the opposite direction, which is almost 
symmetric to the conduction subband about Ep. The first pair of subbands, with three 
constant energy contours within ± 58 meV and a narrow gap E g ~ 8 meV in between them 
(Fig. 2(b)), is special in that it has never appeared in any other stacking configuration. 
The triple-degenerate states are not suitable for low-energy expansion, indicating that the 
effective-mass model can not further deal with the magnetic quantization of these energy 


bands. The second pair of subbands, S^, has a sombrero-shaped and a local energy 
minimum (maximum) and maximum (minimum), situated at around 0.24 (—0.24) eV and 
0.26 (—0.26) eV, respectively. The energy difference between the two extreme points is quite 
narrow, being only about 20 meV. Located away from Ep, the third pair of subbands, S£’ v , 
consists of monotonic parabolic bands with a minimum (maximum) value of about 0.49 
(—0.49) eV. The above-mentioned features of the low-lying energy bands are consistent 
with those by the first-principle calculations (dashed curves), clearly indicating that the 
complex interlayer amomic interactions used in the generalized tight-binding model are 
suitable. 

As for the low-lying band structure, there are important differences among trilayer 
graphenes with distinct stacking configurations. The AAA-, ABA-, and ABC-stacked tri¬ 
layer graphenes possess special band structures with three pairs of linear bands, a pair of lin¬ 
ear bands and two pairs of monotonic parabolic bands, and a pair of sombrero-shaped bands 
and two pairs of parabolic subbands, respectively. The band structures of the AAA and 
ABA stackings can be, respectively, regarded as the superposition of monolayer graphene, 
and monolayer and bilayer graphenes. In contrast, ABC- and AAB-stacking configurations 
exhibit special energy band structures, which are not similar to those of monolayer and 
bilayer systems. It should be noted that the energy band structure of AAB-stacking system 
exhibits an extremely small gap, as indicated in Fig. 2(b). That is to say, while the three 
stacking trilayer graphenes are gapless 2D semimetals with a slight overlap between the 
valence and conduction bands, the AAB-stacked system is a narrow-gap semiconductor. 

A perpendicular magnetic held can quantize electronic states into dispersionless LLs 
with high degeneracy. AAB-stacked trilayer graphene presents rich magnetoelectronic prop- 
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erties. For the {k x = 0 ,k y = 0) state, each LL is two-fold degenerate in the absence of 
spin degeneracy; that is, the wavefunctions localized near 2/6 and 5/6 (4/6 and 1/6) are 
identical. The total carrier density, in which each LL can be occupied, is D—AeBo/hc per 
unit area. There exist three groups of LLs, each group can furthermore be divided into 
two sub-groups corresponding to the two localizations 2/6 and 4/6. Each group consists 
of unoccupied conduction and occupied valence LLs. The former and the latter are almost 
symmetric about Ep and present similar behavior. We first discuss the LLs localized near 
2/6 at Bq = 40 T, being characterized by the subenvelope functions on the six sublattices 
shown in Figs. 3-5. The first, second and third groups are, respectively, initiated at 0, 0.24 
and 0.55 eV for the conduction states, and at 0, -0.23 and -0.58 eV for the valence states. 
Near Ep, the first group starts to show up with the first three LLs (blue lines), which are 
very close to each other. The quantum numbers are determined by the A 1 sublattice with 
the dominating subenvelop function (Fig. 3). Starting from the middle LL, n c { v = 0, which 
is located almost right at Ep, the ordering then increases for the higher conduction states, 
and the lower valence states. In particular, the quantum number of the next conduction 
(valence) LL, which is placed at about E c =10 meV (E v = —5.2 meV), is assigned n\ — 1 
(ril = 1), and so on. For the first group, the quantum number ordering is similar to that 
of monolayer graphene, while the LL spacing is irregular. The largest spacing between LLs 
is that of n c { v = 1 and n c { v = 2, which is about AE C,V = 40 meV; it then decreases with 
increasing quantum numbers. Apparently, there is no simple relationship between the LL 
energy and quantum number. This is in contrast to the relationships found in monolayer 
graphene and 2D electron gas, which can be described by E c,v oc \Jn c ' v and E(n) oc n c,v , 
respectively. 


10 



The destruction of the inversion symmetry in AAB-stacked trilayer graphene results in 
certain important differences for the LL wavefunctions near 4/6 and 2/6, including differ¬ 
ences in the LL spacing, quantum number and spatial distribution of the wavefunctions. 
At B 0 = 40 T, the first, second and third LL groups near the 4/6 center are, initiated at 
0, (0.26 eV, —0.24 eV) and ±0.52 eV, as indicated in Fig. 4(a)-4(c), respectively. For the 
low-lying LLs, it is not easy to define the quantum numbers due to the non-well-behaved 
spatial distributions of the LL wavefunctions. Particularly, the subenvelope functions os¬ 
cillate abnormally as a result of the complex LL anticrossings (details in Fig. 6(c)). The 
eleven LLs nearest to Ep, except for the n = 0 one, are arranged in pairs with very small 
spacings, in which the LL crossings or anticrossings are clearly revealed (Fig. 4(a)). At 
low-lying energy levels, the LL wavefunctions are well-behaved only when the applied field 
strength is sufficiently large (the available magnetic field in experiment is nowadays up to 
80 T [85]), as shown in Fig. 5(a)-5(c) for Bo = 80 T. The dominating subenvelope functions 
of the sublattice B 3 are used to define the quantum numbers. Similarly to those near 2/6 
center, the conduction and valence LLs are also almost symmetric about the n = 0 LL, 
which is located right at Ep. ffowever, the quantum number ordering is slightly differ¬ 
ent to that of monolayer graphene. Particularly, the next unoccupied (occupied) LLs are 
assigned, respectively, n\ = 2 (n\ = 2), n( = 3 n] = 3, = 1 = 1), and so on. 

Concerning the second group, all the LL wavefunctions are single-mode at B 0 = 40 T. 
Also, the LL quantum numbers are determined by the dominating B 1 sublattice for the 
wavefuntions distributed around 2/6 and 4/6 (green color in Figs. 3(b) and 4(b)). The 
assigned quantum numbers of the former are rhf = 1, n^ 1 = 0, n^' = 2, n^ 1 = 3, ... in 
the order of increasing (decreasing) energy for the conduction (valence) LLs, similarly to 
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those of the latter {n^ 0 = 1, n? 2 ' v = 0, n\ 2 V = 2, n 2 V = 3). It should be noted that when 
the held strength is sufficiently large, the quantum number ordering becomes equivalent to 
that of monolayer graphene, as shown in Fig. 5 (b) near the 4/6 center for Bq = 80 T (the 
available magnetic held in experiment is nowadays up to 80 T [85]). On the other hand, 
the third group exhibits a normal quantum number sequence for two localization centers, 
similarly to that of monolayer graphene. That is to say, the ordering is n c ^ v = 0, n 2 v = 1, 
77.3 v = 2, nl’ v = 3, and so on. The LL wavefuntions based on the dominating subenvelope 
function A 2 , are almost identical to those of monolayer graphene, and independent of the 
held strength. 

The F?o-dependent energy spectra are very useful in understanding the rich magnetic 
quantization in AAB-trilayer graphene. The LL energies exhibit monotonic and entangled 
spectra, in which the multi-crossings and -anticrossings occur frequently (Fig. 6). For all 
the LLs of the hrst group localized at 2/6, the quantum numbers have a normal sequence 
when the held strength is very large (B 0 > 40 T). Moreover, the conduction (valence) 
LL energies monotonously decline (grow) when the held strength is reduced to B 0 = 40 
T. With a further decrease of held strength, the L?o-dependent energy spectrum becomes 
oscillatory, which leads to the frequent and pronounced multi-crossings and -anticrossings. 
If two different multi-mode LLs simultaneously possess the same mode, they are forbidden 
to cross each other. For example, the three conduction LLs, nf — 1, n\ — 4 and n\ — 7 
continuously anticross, as indicated in Fig. 6(b). I 11 particular, the anticrossing between 
the former two occurs in the range of Bq ~ 19 — 21 T, while in the latter two they occur 
between B 0 ~ 9 — 11 T; the anticrossing between the hrst and the third conduction LLs 
appears at Bo ~ 11 — 13 T. The LL anticrossings clearly indicate the unusual sequence of 


12 


quantum numbers. For B 0 > 21 T, the main n\ = 1 mode has a normal distribution on 
the dominating sublattice, in which the side modes of n\ = 4 and n\ = 7 are faint. With 
the decreasing of Bo, the nf — 1 main mode declines and the nl = 4 side mode quickly 
reaches the maximum at B 0 = 20 T in the range of E c ~ 20 — 40 meV. At the center 
of the anticrossing regions, their subenvelop functions have the same n — 1 and n = 4 
modes, which forbids them to have the same energy. Similarly, the n\ = 4 and n\ = 7 LLs 
anticross in the range of E c ~ 15 — 35 meV. The anticrossing center is at Bq ~ 10 T where 
the comparable n = 4 and n = 7 modes exist in two subenvelope functions. Moreover, the 
other conduction LLs exhibit similar f? 0 -dependent energy spectra. Particularly, there are 
other pairs of LLs, such as n\ = 2 and n\ = 5, n\ = 3 and n'f = 6 , ... , for which the LL 
anticrossings appear, respectively, in the ranges of 16 T < B 0 < 18 T ( E c ~ 18 — 38 meV), 
13 T < B 0 < 15 T ( E c ~ 17 — 37 meV), and so on. The LL anticrossings form a wide 
stateless region which is indicated in Fig. 6 (b). Generally, the main modes of anticrossing 
LLs are different by 3 at large field strengths and 6 at small ones. 

Furthermore, there exists another region of LL anticrossings which is indicated in Fig. 
7(a). These anticrossings happen at Bo < 19 T and form very narrow in-between gaps. 
For example, the n\ — 1 and n\ = 7 LLs are forbidden to cross in the energy range of 
E c ~ 34 — 36 meV (Fig. 7(a)). Likewise, the other pairs of LLs, such as (n'f = 2 and 
n\ = 8 ) ( E c ~ 39 — 40 meV), {n\ = 3 and nl = 9 ) ( E c ~ 41 — 41.5 meV), anticross at 
about B 0 = 11 T and B 0 = 10.6 T, respectively. With a further decrease of L? 0 , the = 4 
and nl = 10 ( n\ = 5 and n'f = 11) LLs anticross each other at around 5o = 6T (Bo = 5 T) 
in the range of E c ~ 36 — 38 meV ( E c ~ 40 — 42 meV), and so on. In short, for each pair 
of anticrossing LLs, the main modes differ by 6. In addition to the intragroup anticrossings 
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of the LLs in the first group, the second group LLs also avoid crossing each other at smaller 
field strengths (Bq < 5 T) due to the sombrero-shaped energy dispersions [6}[38] shown in 
Fig. 7(b). The third group consists of monotonic LLs without intragroup anticrossings. 

The intergroup LL anticrossings occur among all three groups at sufficiently large field 
strengths. For instance, the n\ = 3 LL anticrosses the = 0 and = 3 ones continuously 
at B 0 ~ 75 — 97 T, in the range of, respectively, E v ~-23—18 meV and E v ~-27—23 
meV, as shown in Fig. 7(c). In general, the main modes of anticrossing LLs are different 
by 3m, with m being an integer. It should be noted that the interlayer atomic interaction 
y5 between the non-vertical sites in the second and third layers induces both intergroup 
and intragroup LL anticrossings. This is similar to ABC-stacked trilayer graphene [ 38] . 

Significantly, near E Fl the n = 0 LL presents the unusual field-dependent LL wavefunc- 
tions, as indicated in Fig. 6(b). When the magnetic field is sufficiently large, e.g., B 0 > 22 
T, both n — 0 and n c,v = 3 LLs present the well-behaved wavefunctions. However, the 
n = 0 LL is forbidden to cross the n c = 3 one when the field strength is reduced to B 0 ~ 22 
T. Accordingly, there appear some n = 0 mode wavefunctions near the turning point of 
the n c = 3 LL, which is located at Bq ~ 20 T and E c ~ 13 meV. Besides, the n c = 3 
and n v = 3 LLs also weakly anticross right at their extreme points, forming an exceedingly 
narrow gap in the anticrossing center. As the field strength is gradually decreased, the 
n = 0 and n c = 3 LLs anticross each other again at B 0 ~ 16 T. Therefore, three LLs in 
the range of E c ~-10—15 meV and B 0 ~ 16 — 22 T exhibit an entangled energy spectrum 
including the special anticrossings. 

On the other hand, there exist certain critical differences between the field-dependent 
LL spectra at the 4/6 and 2/6 localization centers. For the first group, the intragroup 
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LL anticrossings of the former appear in a wider magnetic-field range, 0 T < B 0 < 80 T 
(Fig. 6(c)). Specifically, the complex LL anticrossings of the n — 0 and n c,v = 3 LLs are 
revealed in the range of larger held strengths (Bo ~ 32 — 45 T) compared to that of the 
latter (B 0 ~ 16 — 22 T), as indicated in Figs. 6(b)-(c). The n = 0 LL anticrosses the 
n v = 3 LL twice, leading to the fact that the zero-mode LLs only exist above Ep, which 
is in contrast to those near the 2/6 center. Moreover, the stateless gap in the anticrossing 
center of n c = 3 and n v — 3 LLs is about 5 meV, clearly larger than that at 2/6 localization. 

The I?o-dependent LL wavefunctions of AAB-stacked trilayer graphene exhibit diverse 
real-space distributions. They can be verified by the spectroscopic-imaging STM measure¬ 
ment, as done for 2D electron gas [83] and topological insulators [84]. From the mea¬ 
surements on the variations of the local density of states (DOS) in graphene planes, STM 
reveals a nodal structure corresponding to the well-behaved LL wavefunctions. For the 
hybridized wavefunctions due to the LL anticrossings in AAB-trilayer graphene, the results 
from the STM measurements are expected to be helpful in distinguishing the main mode 
and the side modes of the LLs. Moreover, the energy spectra and the internal structure of 
the wavefunctions could be directly examined by STM and STS measurements (discussed 
later), respectively, giving a useful identifiable picture of the LLs. 

Reflecting the main characteristics of LL spectra, the DOS is defined as D(uj ) = 
Y h s tBZ — k))dk. The discrete LLs lead to many symmetric delta-function-like 

n c ,n v 

peaks in the DOS, where the peak intensities are proportional to the LL degeneracy [861(87] . 
Such peaks correspond to the sharp structures in the differential conductance map of 
dl/dV—V from STS measurements. For the DOS of the AAB-stacked trilayer graphene at 
B=40 T, the first, second and third groups, respectively, have the onset conduction-states 
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about 0, 0.26 and 0.52 eV, as shown in Fig. 8(a). Each group consists of two sub-groups, 
which are separated by a symmetry-breaking induced energy difference of about 10 meV. In 
the vicinity of the Fermi level, the peaks of the first groups are neither regularly sequenced 
nor do they follow a simple relationship as for monolayer graphene. This is due to the 
fact that the special quantization effects on the lowest subband leads to the multi-crossings 
and anticrossings within the band width. On the other hand, the second and third groups 
exhibit a normal sequence in the order of increasing energies. A crossover of LLs in different 
groups induces higher DOS, which is thus expected to exhibit stronger tunneling currents 
in experimental measurements. 

The DOS exhibits distinct characteristics among different stacking configurations, such 
as the peak sequence, intensity, energy and splitting. In AB-stacked trilayer graphene, 
peaks that follow the sequence in monolayer and bilayer graphenes are observed, as pre¬ 
sented in Fig. 8(b). Similarly, for the AAB stacking, the half-intensity peaks near the 
Fermi level are attributed to the symmetry-broken structure, which leads to a considerable 
energy splitting of about 10 meV. However, this splitting is hardly observable at higher en¬ 
ergies. The AA-stacked trilayer graphene exhibits three groups of monolayer-like sequence 
of peaks [851ISS1IS5] located at energies described by the simple relationship E c,v oc \/n c ’ v B, 
as indicated in Fig. 9(a). On the contrary, the DOS of ABC-stacked trilayer graphene is 
distinct from that of monolayer and bilayer graphenes, as shown in Fig. 9(b). The excep¬ 
tionally high peak at the Fermi level is a superposition of three peaks corresponding to the 
Dirac points; its intensity is proportional to the number of graphene layers. The essential 
differences of the DOS can be verifed by STS [86TT89] ; those profiles can then be used as a 
tool to identify the stacking configuration of graphene sheets. 
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4. Conclusion 


The generalized tight-binding model is developed to investigate the rich magneto- 
electronic properties of AAB-stacked trilayer graphene. This system exhibits special energy 
bands and thus a rich magnetic quantization. The complex interlayer atomic interactions 
indicate that the low-energy expansion about the K point is not suitable in obtaining energy 
dispersions; therefore, the effective-mass model can not be used to achieve further magnetic 
quantization. The three zero-field pairs of energy bands contain oscillatory, sombrero¬ 
shaped and parabolic dispersions, which could be examined by ARPES [82]. They are 
magnetically quantized into three groups of LLs, defined by the dominating subenveloped 
functions. The field-dependent energy spectra, in which the frequent intragroup and inter¬ 
group LL anticrossings happen simultaneously, are very complex. There exist important 
differences between the LLs near 4/6 and 2/6 in terms of the LL splitting, quantum num¬ 
ber, spatial distribution of wavefunctions, and anticrossing LL energy spectra. Moreover, 
the electronic properties of AAB-stacked trilayer graphene sharply contrast to those of the 
AAA-, ABA- and ABC-stacked ones, such as state degeneracy, LL anticrossing behavior, 
initial energies of each group, quantum number ordering and wavefunction distributions. 
This principally demonstrates the unique characteristics of the special AAB configuration. 

The first and the second LL groups exhibit unusual energy spectra with continuous 
intragroup LL anticrossings. Especially for the first group, each low-lying LL is forbidden 
to cross with the others more than once. The LL anticrossing region near Ep has a wide 
stateless gap, while that of the other region is quite narrow. The main modes of LLs during 
anticrossings are different by 3 and 6 for the former, and only 6 for the latter. Furthermore, 
the intergroup LL anticrossings among three groups are revealed only at sufficiently large 
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field strengths. Remarkably, the DOS exhibits many symmetric delta-function-like peaks 
with an irregular sequence near Ep due to the LL multi-anticrossings. The web-like energy 
spectra and the seriously distorted wavefunctions could be examined by STS and STM 
experimental measurements, respectively. The complex LL spectra are expected to induce 
the rich mechanical, excitonic, thermal and optical properties of the AAB-stacking system. 
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FIGURE CAPTIONS 


Fig. 1 - The interlayer atomic interactions and the geometric structure under a uniform 
magnetic field Bqz. The shaded region corresponds to a rectangular unit cell. The first 
and second layers have the same (x, y) projections. 

Fig. 2 - The (a) energy band structures of AAB-stacked trilayer graphene with (b) a 
narrow energy gap. 

Fig. 3 - The distributions of the LL wavefunctions based on the dominating sublattices 
centered at the (a)-(c) 2/6 localization under Bo = 40 T. 

Fig. 4 - The distributions of the LL wavefunctions based on the dominating sublattices 
centered at the (a)-(c) 4/6 localization under B 0 = 40 T. 

Fig. 5 - The distributions of the LL wavefunctions based on the dominating sublattices 
centered at the (a)-(c) 4/6 localization under B 0 = 40 T. 

Fig. 6 - Three groups of field-dependent LL energy spectrum for (a) the 2/6 localization 
center; the low-lying LL spectra for (b) the 2/6 and (c) 4/6 centers. 

Fig. 7 - The (a) second region of intragroup LL anticrossings in the first group, the (b) 
frequent intragroup LL anticrossings of the second group at low field strength, and the (c) 
intergroup LL anticrossings of the first and second groups. 

Fig. 8 - The DOS of (a) AAB- and (b) ABA-stacked trilayer graphene at B 0 = 40 T. 

Fig. 9 - The DOS of (a) AAA- and (b) ABC-stacked trilayer graphene at B 0 = 40 T. 
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Figure 1: The interlayer atomic interactions and the geometric structure under a uniform 
magnetic held B 0 z. The shaded region corresponds to a rectangular unit cell. The first 
and second layers have the same (x, y) projections. 
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Figure 2: The (a) energy band structures of AAB-stacked trilayer graphene with (b) a 
narrow energy gap. 
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Figure 3: The distributions of the LL wavefunctions based on the dominating sublattices 
centered at the (a)-(c) 2/6 localization under B 0 = 40 T. 
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Figure 4: The distributions of the LL wavefunctions based on the dominating sublattices 
centered at the (a)-(c) 4/6 localization under B 0 = 40 T. 
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Figure 5: The distributions of the LL wavefunctions based on the dominating sublattices 
centered at the (a)-(c) 4/6 localization under B 0 = 80 T. 
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Figure 6: Three groups of field-dependent LL energy spectrum for (a) the 2/6 localization 
center; the low-lying LL spectra for (b) the 2/6 and (c) 4/6 centers. 
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Figure 7: The (a) second region of intragroup LL anticrossings in the first group, the (b) 
frequent intragroup LL anticrossings of the second group at low field strength, and the (c) 
intergroup LL anticrossings of the first and second groups. 
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Figure 8: The DOS of (a) AAB- and (b) ABA-stacked trilayer graphene at B 0 = 40 T. 
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Figure 9: The DOS of (a) AAA- and (b) ABC-stacked trilayer graphene at B 0 - 


40 T. 
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